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Abstract 

A family of models of counterrotating and rotating relativistic thin disks of infinite exten- 
ts , sion based on a charged and magnetized Kerr-NUT metric are constructed using the well-known 
"^ I "displace, cut and reflect" method extended to solutions of vacuum Einstein-Maxwell equations. 
• ■ The metric considered has as limiting cases a charged and magnetized Taub-NUT solution and 
L». , the well known Kerr-Newman solutions. We show that for Kerr-Newman fields the eigenval- 
00 I ues of the energy-momentum tensor of the disk are for all the values of the parameters real 
^^ ' quantities so that these disks do not present heat flow in any case, whereas for charged and 
magnetized Kerr-NUT and Taub-NUT fields we find always regions with heat flow. We also 
>v>( I find a general constraint over the counterrotating tangential velocities needed to cast the surface 
^H ' energy-momentum tensor of the disk as the superposition of two counterrotating charged dust 
fluids. We show that, in general, it is not possible to take the two counterrotating fluids as 
circulating along electrogeodesics nor take the two counterrotating tangential velocities as equal 
and opposite. 
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1 Introduction 

Stationary or static axially symmetric exact solutions to the Einstein field equations representing rela- 
tivistic thin disks are of great astrophysical importance since they can be used as models for certain 
galaxies, accretion disks, and the superposition of a black holes and a galaxy or an accretion disk as 
in the case of quasars. Disk sources for stationary axially symmetric spacetimes with electromagnetic 
fields, especially magnetic fields, are also of astrophysical importance in the study of neutron stars, 
white dwarfs and galaxy formation. In such situation one has to study the coupled Einstein-Maxwell 
equations. 

Exact solutions which describe relativistic static thin disks without radial pressure were first 
studied by Bonnor and Sackfield [1], and Morgan and Morgan [2j, and with radial pressure by 
Morgan and Morgan [3]. Also thin disks with radial tension were considered [1]. Several classes 
of exact solutions of the Einstein field equations corresponding to static thin disks with or without 
radial pressure have been obtained by different authors [5]-[T2]. Rotating thin disks that can be 
considered as a source of a Kerr metric were presented by Bicak and Ledvinka [13], while rotating 
disks with heat flow were were studied by Gonzalez and Letelier [H]. The nonlinear superposition 
of a disk and a black hole was first obtained by Lemos and Letelier [9]. Perfect fluid disks with 
halos were studied by Vogt and Letelier p^. The stability of some general relativistic thin disks 
models using a first order perturbation of the energy-momentum tensor was investigated by Ujevic 
and Letelier [TB]. 

Thin disks in precense of electromagnetic field have been discussed as sources for Kerr- Newman 
fields [TTlllH], magnetostatic axisymmetric fields [19], conformastationary metrics [20], while models 
of electrovacuum static counterrotating dust disks were presented in pT] . Charged perfect fluid disks 
were also studied by Vogt and Leteher [22] , and charged perfect fluid disks as sources of static and 
Taub-NUT-type spacetimes by Garcia- Reyes and Gonzalez [23l El] . 

In all the above cases, the disks are obtained by an "inverse problem" approach, called by Synge 
the "5^-met/iocf ' [25] . The method works as follows: a solution of the vacuum Einstein equations is 
taken, such that there is a discontinuity in the derivatives of the metric tensor on the plane of the disk, 
and the energy-momentum tensor is obtained from the Einstein equations. The physical properties 
of the matter distribution are then studied by an analysis of the surface energy-momentum tensor 
so obtained. Another approach to generate disks is by solving the Einstein equations given a source 
(energy-momentum tensor). Essentially, they are obtained by solving a Riemann-Hilbert problem 
and are highly nontrivial [26]-[32]. A review of this kind of disks solutions to the Einstein-Maxwell 
equations was presented by Klein in [33] . 

Now, when the inverse problem approach is used for static spacetimes, the energy-momentum 
tensor is diagonal and its analysis is direct and, except for dust disks, the solutions obtained have 
anisotropic sources with azimuthal stress different from the radial stress. On the other hand, when 
the considered spacetime is stationary, the obtained energy-momentum tensor is non-diagonal and 
the analysis of its physical content is more involved and, in general, the obtained source is not only 
anisotropic but with nonzero heat flow. Due to this fact, there are very few works about of stationary 
disks and they are limited to disks without heat flow [131 [IZl [21] . Indeed, only in one work disks with 



heat flow has been considered, but only a partial analysis of the corresponding energy-momentum 
tensor was made [H]. 

The above disks can also be interpreted as made of two counterrotating streams of freely moving 
(charged) particles, i.e., which move along (electro-)geodesics, as was also indicated in p/7j. This 
interpretation is obtained by means of the Counterrotating Model (CRM) in which the energy- 
momentum tensor of the source is expressed as the superposition of two counterrotating fluids. Now, 
even though this interpretation can be seen as merely theoretical, there are observational evidence 
of disks made of counterrotating matter as in the case of certain 5*0 and Sa galaxies [35]-[38]. Has 
also been observed in some spiral galaxies as in NGC3626 the presence of counterrotating ionized 
(charged) gas disks [39]. Indeed recent investigations have shown that there is a large number of 
galaxies (the first was NGC4550 in Virgo) [38j which show counterrotating streams in the disk with 
up to 50% counterrotation. It is believe that the presence of counterrotating matter components in 
these galaxies is the consecuence of the accretion or merger of galaxies. Thus these counterrotating 
streams may be the result from the capture by a massive early-type galaxy of a gas-rich dwarf galaxy 
that was orbiting in the opposite sense to the rotation of the main galaxy. 

The purpose of the present paper is twofold. In first instance, we present a analysis of the energy- 
momentum tensor and the surface current density for electrovacuum stationary axially symmetric 
relativistic thin disks of infinite extension without radial stress and when there is heat flow. And, 
in the second place, we present a study of the counterrotating model (CRM) for these stationary 
thin disks. The paper is structured as follows. In Sec. II we present a summary of the procedure to 
obtain models of rotating thin disks with a purely azimuthal pressure and currents, using the well 
known "displace, cut and reflect" method extended to solutions of Einstein-Maxwell equations. In 
particular, we obtain expressions for the surface energy-momentum tensor and the surface current 
density of the disks. 

In Sec. Ill the disks are interpreted in terms of the CRM. We find a general constraint over 
the counterrotating tangential velocities needed to cast the surface energy-momentum tensor of the 
disk as the superposition of two counterrotating charged fluids made of dust or pressureless matter 
(collisionless particles). We show that, in general, it is not possible to take the two counterrotating 
fluids as circulating along electro-geodesies nor take the two counterrotating tangential velocities as 
equal and opposite. 

In the following section. Sec. IV, a family of models of counterrotaing and rotating relativistic 
thin disks based on a magnetized and charged Kerr-NUT metric is considered. This metric has 
as limiting cases a charged and magnetized Taub-NUT solution and the well known Kerr-Newman 
solutions. We show that for Kerr-Newman fields the eigenvalues of the energy-momentum tensor of 
the disk are for all the values of the parameters real quantities so that this disks do not present heat 
flow in any case, whereas for charged and magnetized Kerr-NUT and Taub-NUT fields we find always 
regions with heat flow. We then analyze the CRM for these disks and we study the angular velocities, 
surface energy densities and electric charge densities of both streams when the two fluids move along 
electrogeodesics and when they move with equal and opposite velocities. Also the stability against 
radial perturbation is analyzed in all the cases. Finally, in Sec. V, we summarize our main results. 
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2 Electrovacuum rotating relativistic thin disks 

A sufficiently general metric for our purposes can be written as the Weyl-Lewis-Papapetrou line 
gIgixigiiIj i341 

where we use for the coordinates the notation {x^ , x^ , x"^ , x^) = {t,ip,r,z), and \E', W, and A are 
functions of r and z only. The vacuum Einstein-Maxwell equations, in geometric units in which 
SnG = c = fiQ = Eq = 1, are given by 

(2a) 
(2b) 

with 

(3a) 
(3b) 

where Aa = {At, A^, 0, 0) and the electromagnetic potentials At and A^ are also functions of r and 
z only. 

For the metric ([T]), the Einstein-Maxwell equations are equivalent to the system [40j 

V ■ [r-V(VA^ - W VAt] = 0, (4a) 

V ■ [f-'VAt + r-^fW{VA^ - W VA] = 0, (4b) 

V ■ [r-^fVW - 2r-^fAt{VA^ - W VAt)] = 0, (4c) 
/VV = V/ ■ V/ - r-^fVW ■ VW + fWAt ■ VAt 

+ r-^f{VA^ - WVAt) ■ {VA^ - WVAt), (4d) 

A. = r(v&^.-n)-^(W,i-VV,^Je^*-l(rV2^->VV^)(4.-4J 

+ 1«, - AlJe'" - \yV{A^,rAt,r - A^,.A.)e'^ (4e) 

A,, = 2rM/,,M/,,-;iw,.W..e^*-i(rV2'^-WV*)A,,A., 

+ ^A^^rA^^.e^"" - ^WiA^,rAt,, + A^^,At^r)e^^ , (4f) 

where V is the standard differential operator in cylindrical coordinates and / = e^*. 

In order to obtain a solution of (!2a|) - (l2bl) representing a thin disk at 2; = 0, we assume that the 
components of the metric tensor are continuous across the disk, but with first derivatives discontin- 
uous on the plane 2 = 0, with discontinuity functions 



Thus, by using the distributional approach [HI |l2l |33] or the junction conditions on the extrinsic 
curvature of thin shells [HlHSlllGj, the Einstein-Maxwell equations yield an energy- momentum tensor 
Tab = Tf^ + TX\ where T^' = Qab S{z), and a current density Ja = jaS{z) = -2e^^''-''^Aa,J{z), 
where S{z) is the usual Dirac function with support on the disk. T^^™ is the electromagnetic tensor 
defined in Eq. ( I3al) . ja is the current density on the plane z = 0, and 

is the distributional energy-momentum tensor. The "true" surface energy-momentum tensor (SEMT) 
of the disk, Sah, and the "true" surface current density, j^, can be obtained through the relations 



Sab = j TT dsn = e^-* Qab , (6a) 

ia = Ja dSn = 6^'^ ja , (6b) 



where dsn = y^gll dz is the "physical measure" of length in the direction normal to the disk. 
For the metric ([1]), the nonzero components of 5*^ are 



g^-A 



Si = '—^ [2r2(A„ - 2V1/,, ) - e^^WW,, ] , (7a) 

5? = -£!^[4r2>Vvl>„+(r2+ >VV*)W,.], (7b) 

^0^ = ^[e^^W,.], (7c) 

S\ = ^[2r2A„+e^*WW,.], (7d) 
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and the nonzero components of the surface current density j^ are 

j, = -2e*-^A,., (8a) 

U = -2e*-M^,„ (8b) 

where all the quantities are evaluated at 2; = 0^. 

These disks are essentially of infinite extension. Finite disks can be obtained introducing oblate 
spheroidal coordinates, which are naturally adapted to a disk source, and imposing appropriate 
boundary conditions. These solutions, in the vacuum and static case, correspond to the Morgan 
and Morgan solutions [2]. A more general class of solutions representating finite thin disks can be 
constructed using a method based on the use of conformal transformations and solving a boundary- 
value problem [3l El H [121 [23l [2l] . 

Now, in order to analyze the matter content of the disks is necessary to compute the eigenvalues 
and eigenvectors of the energy-momentum tensor. The eigenvalue problem for the SEMT (I7a|) - (]7dp 



s^ e =x r, (9) 



has the solutions 

A± = 1(t±v^), (10) 

where 

T = S', + Si, D= (Si - S',y + 4 ^0 S'„ (11) 

and Xr = Xz = 0. For the metric (IT]) 



;2{*-A) 

T = 4e*-^(A,,-vI/,,), (13) 



D = A^-^-iAr'^l-Wy) = A''-B^, (12) 



where 

A = 4^ ,e*-^, S = -W.e^*-^. (14) 

r 

The corresponding eigenvectors are 

a = ( e^, e±, 0, 0), 

^a _ e*-^(0, 0,1,0), (15) 

ya ^ e*-'^(0, 0,0,1), 

with 

to 



,K,.«J^2A^.e-^^^-^^J. (16) 

where 

N± = \nD{-VD±A). (17) 

We only consider the case when D > 0, so that the two eigenvalues X± are real and different 
and the two eigenvectors are orthogonal, in such a way that one of them is timelike and the other is 
spacelike. Since |y4| > a/D, from (IT71) follows that when A > the negative sign corresponds to the 
timelike eigenvector while the positive sign to the spacelike eigenvector. When A < we have the 
opposite case. So the function \1/ ^ determines the sign of the norm. 

Let V"" be the timelike eigenvector, VaV"" = —1, and W"" the spacelike eigenvector, WaW"" = 1. 
In terms of the orthonormal tetrad or comoving observer ea'' = {V"^, W'', X^, Y^}, the SEMT and the 
surface electric current density may be decomposed as 

(18a) 
(18b) 

where 

(19) 



Sab 


= eVaVb+p^WaWh 


Ja 


= fVa+fWa, 


€ = 


-A±, p^ = A^, 



are, respectively, the surface energy density, the azimuthal pressure, and 

f = -VX, j' = WX, (20) 

are the surface electric charge density and the azimuthal current density of the disk measured by 
this observer. In (iT9l) the sign is chosen according to which is the timelike eigenvector and which is 
the spacelike eigenvector. However, in order to satisfy the strong energy condition g = e + p^p > 0, 
where g is the effective Newtonian density, we must choose ^_ as the timelike eigenvector and ^+ 
as the spacelike eigenvector. These condition characterizes a disk made of matter with the usual 
gravitational attractive property. Consequently \E' ^ must be taken positive. So we have 



-A-, p^ = X+, (21) 



and 



ue 



V = -^^{Sl-Sl-^\ (22a) 

= ^'^^ Si (22b) 

where z/ = ±1 so that the sign is chosen according to the causal character of the timelike eigenvector 
(observer's four- velocity), 

iy° = ^2=^° (23a) 

W = -l^{Sl-S', + ^/D), (23b) 



where 



M = Vd lgnVD + 2rWB + {r^e-^"^ + W^e^'^)^} . (24) 



When D < 0, the eigenvalues X± and the eigenvectors ^± are complex conjugates, A_ = A+, 
^a _ ^a_ jj^ ^ comoving orthonormal tetrad e^^ = {V^,W^,X'',Y''} the eigenvectors C,± can be 
expressed as C,± = V"" ± iW"". So we can write the SEMT in the canonical form 

Sab = eVaVt, + q{VaWf, + WaV,)+p^WaWh, (25) 

so that it can be interpreted as the energy-momentum tensor of a matter distribution with propa- 
gation of heat in the tangential direction. The energy density, the azimuthal pressure, and the heat 
flow function are respectively 

^ = - 2 ' ^- = 2 ' ^ = ^— ^^^^ 



3 Counterrotating charged dust disks 

We now consider, based on Refs. p^ and [23j, the possibility that the SEMT S"'^ and the current 
density f can be written as the superposition of two counterrotating charged fluids that circulate in 
opposite directions; that is, we assume 

gab ^ gab ^ gab ^ ^27a) 

/ = J++J-, (27b) 

where the quantities on the right-hand side are, respectively, the SEMT and the current density of 
the prograde and retrograde counterrotating fluids. 

Let f/|. = {Uj., U^, 0, 0) = t/±(l, uj±, 0, 0) be the velocity vectors of the two counterrotating fluids, 
where uj± = U]_/U^ are the angular velocities of each stream. In order to do the decomposition fl27al) 
and (I27bp we project the velocity vectors onto the tetrad ea'', using the relations [17] 



Ui = e\Ui, Ul = ey^Ul. (2 

In terms of the tetrad (fT5|) we can write 



ja 



-± 


vl-4 ' 


VI- 


- v'iv+U't - \/l - v\v_Ul 


f + — f _ 


v/1- 


- vlUl - Vl - vlUl 



(29) 



so that 

ya ^ x^ ■^-■^^-^^ v^ "+"~^ + , (30a) 

w = yr^^T - yr^^-_, (30b) 

f + — v^ 

where v± = U^/U^ are the tangential velocities of the streams with respect to the tetrad. 

Another quantity related with the counterrotating motion is the specific angular momentum of 
a particle rotating at a radius r, defined as h± = g^aU±. This quantity can be used to analyze the 
stability of circular orbits of test particles against radial perturbations. The condition of stability, 

\^ > 0, (31) 

ar 

is an extension of Rayleigh criteria of stability of a fluid in rest in a gravitational field [18|. For an 
analysis of the stability of a rotating fluid taking into account the collective behavior of the particles 
see for example Refs. [HI [16] . 



Substituting fl30ap and (]30bp in (125|1 we obtain 



S' 



ab 



+ 



F{v_,v_){l-vl) U^Ul 
F{v+,v+){l-vl) U^Ut 

(f+ — f_)^ 



where 



F(t;i, f2) = tviV2 - g(t;i + t;2) + P^- (32) 

Clearly, in order to cast the SEMT in the form fl27ap . the mixed term must be absent and therefore 
the counterrotating tangential velocities must satisfy the following constraint 



Fiv+.v. 



0, 



(33) 



where we assume that \v±\ ^ 1. 

Then, assuming a given choice for the tangential velocities in agreement with the above relation, 
we can write the SEMT as f l27ap with 



a rrb 



5f = e± U^U'^ 



(34) 



so that we have two counterrotating dust fluids with energy densities, measured in the coordinates 
frames, given by 



e± 



4 



(e%-g), 



(35) 



.% -v±_ 

Thus the SEMT S""^ can be written as the superposition of two counterrotating dust fluids if, and 
only if, the constraint (133|) admits a solution such that f+ 7^ t'-. 

Similarly, substituting fl30ap and (JSObp in (]18bp we can write the current density as (]27bp with 



J± = (^±U± 



(36) 



where a± are the counterrotating electric charge densities, measured in the coordinates frames, which 
are given by 



cr± 



4 



v± -% 



(/-/%)• 



(37) 



Thus, we have a disk makes of two counterrotating charged dust fluids with energy densities given 
by (135|1 . and electric charge densities given by (!37|) . 

As we can see from Eqs. (!29|) . (!35l) and (1371) . all the main physical quantities associated with the 
CRM depend on the counterrotating tangential velocities v±. However, the constraint (133|1 does not 



determine v± uniquely so that we need to impose some additional requirement in order to obtain a 
complete determination of the tangential velocities leading to a well defined CRM. 

A possibility, commonly assumed [T71I3T]. is to take the two counterrotating streams as circulat- 
ing along electrogeodesics. Now, if the electrogeodesic equation admits solutions corresponding to 
circular orbits, we can write this equation as 



le±gab,rU^Ui = -a^FrJJl. (38) 



In terms of uj± we obtain 



2^±(f^±)^(fl'll,rt^| + '^g01,r^± + fi'OO.r) = -<y±U±{A,r + A^,r^±)- (39) 



From (l27aD . (l27bD . ([21, and ([2SD we have 



a^Ul = L_^^ (40a) 

COI _ , , cOO 

e±{Ulf = ^ ^^i^, (40b) 

and substituting (I40ap and fl40bp in fl39|) we find 

-(/ - ^T/)(A,r + A^,r-W±), (41) 

and using (I40cj) we obtain 

-[5°Y - S''f + co^iS^'f - 5°V)](A,. + A^,.^±). (42) 

Therefore we conclude that 

J-l 



with 



,■0^01 _ ji5'00 
7^1 = gn,r + 2A^,r ^01^01 _ 500511 ' (44^) 

jOgoi_jigoo /goi-jOgii 

-'2 5'01,r + ^t.r- ^015-01 _ 5-005-11 "^ "^'/'.r- ^01 501 _ 500511' '^ "'' 

/I5OI _ 70511 

Ts = ^oo.r + 2A,. go,go,_googi, . (44c) 

10 



It is easy to see when D > that electrogeodesic velocities satisfies the constraint ( 133|) . In the 
fact, in terms of u!± we get 

~Wo + Wiuj± 



v± = - 



and, by using fH3|) . we have that 



v^v_ 



Vo + Viu± 
TiWj - 2T2W0W1 + Tgiyf 



(45) 



(46) 



so that, using (llSal) . we get 



F{v+,v. 



32e4{*-A)^2 /^2^ ,x/D + Ar^A ,^ .e*"^ - W^e^*- 



r^{A + VD)p^{S^'S^' - S^^S'^){T^V,^ - 2T2V0V, + T^V,^) 



2r T 



Y^'e''')A,,A 



t,z 



A^^rA^^z^ H yyyA^^rAt^z + A^^zAt^r)^ 



2* 



(47) 



Finally, using the Einstein-Maxwell equation f l4fl) follows immediately that F(f+,w_) vanishes and 
therefore the electrogeodesic velocities satisfy the constraint ( l33l) and so, if the electrogeodesic equa- 
tion admits solutions corresponding to circular orbits, we have a well defined CRM. When there 
is heat fiow a proof as the previous one is not trivial and in this case we take the counterrotating 
hypothesis (!34ll and (J36ll as an ansatz. 

Another possibility is to take the two counterrotating fiuids not circulating along electrogeodesics 
but with equal and opposite tangential velocities. 



v± = ±v = ±\ p^/ 



(48) 



This choice, that imply the existence of additional interactions between the two streams (e.g. col- 
lisions), leads to a complete determination of the velocity vectors. However, this can be made only 
when < p<^/e < 1. So in precense of heat fiow by Eq. fl2Bl) we have p<^/e = —1 which corresponds 
to an imaginary velocity and therefore these disks have an unphysical behavior in the regions where 
there is heat fiow. 

In the general case, the two counterrotating streams circulate with different velocities and we can 
write ([33]) as 



V+V- = -. 

e 



(49) 



However, this relation does not determine completely the tangential velocities, and therefore the 
CRM is undetermined. In summary, the counterrotating tangential velocities can be explicitly de- 
termined only if we assume some additional relationship between them, like the equal and opposite 
condition or the electro-geodesic condition. Now, can happen that the obtained solutions do not 
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satisfy any of these two conditions. That is, the counterrotating velocities are, in general, not com- 
pletely determined by the constraint (133|) . Thus, the CRM is in general undetermined since the 
counterrotaing energy densities and pressures can not be explicitly written without a knowledge of 
the counterrotating tangential velocities. 



4 Disks from a Charged and magnetized Kerr-NUT solution 

As an example of the above presented formalism, we consider thin disk models obtained by means of 
the "displace, cut and reflect" method applied to the electromagnetic generalization of the Kerr-NUT 
metric, which can be written as 



A 
W 

At 



-In 
2 

2 

2kc 



a^x^ + hiy^ — c 



V? + f ^ 

,2„.2 



a{x^ + h{y^ — c 
hi{l-y^) 



afx^ 



v^ 



+ 6fy2 _ c2) 

a2U + h2V 
V? + v"^ 



\a1a2X 



hih2y + -c{l + c^)]+h2y 



, v/2(c2 - 1) r {a2U + h2v)[-y{hiy + 2062) + h] , , 

-k- <; — + 622/ 



Qi 



u^ 



where 



u 



aix + ca2. 



biy + cb2 



(50a) 
(50b) 
(50c) 
(50d) 
(50e) 

(51) 



and af + bf = a2 + b2 = c^ > 1, being c the parameter that controls the electromagmetic field, x and 
y are the prolate spheroidal coordinates which are related with the Weyl coordinates by 

r^ = k'^{x^ -l){l-y^), z + zo = kxy, (52) 

where l<a;<oo,0<y<l, and k is an arbitrary constant. Note that we have displaced the 
origin of the z axis in zq. This solution can be generated, in these coordinates, using the well-known 
complex potential formalism proposed by Ernst |10] using as seed solution the Kerr-NUT vacuum 
solution [33]. So when c = 1 this solution reduces to the Kerr-NUT vacuum solution. When 61 = 
we have a charged and magnetized Taub-NUT solution [23] , and when 62 = we have the well known 



12 



Kerr-Newman solution. Let T = kT, D = k^D, j^ = kj^ be, therefore 

^ = (3.2 _ ^2)3/2^^2 + ^2)3/2 {-QlQ2p' + [-2cf - Sb^b^f 

+c(l — c^)y + bib2]x^ + 3aia2y{l — y^)x^ 

+ [-bib2f + c(l - c^)f + 36162^ + 2c^]xy + aia2f] , (53a) 

^ = T^ ^:r--^—^r:T{bj{c^x'^-al)f + 2bib2ic^x + cx + 2aia2)xf 

[x^ — y^jyu^ + v'^j'^ ^ 

+ [c^a\x^ + 2aia2c{(? + l)(a;^ + l)a; — (3c^a^ + ?)(?a\ — Qa\a\ — c® 

-Ac^ - c^)x^ + c^aa]^^ + 2bib2{2aia2X + c^ + c)x^y - bl{alx^ - c^)x^] , (53b) 

~h = - (^2 _ ^2)1/2(^2 /^2)5/2 {blhxf + &^(3aia2X^ + 2c'x - a^a2)f 

—bib2{Sa1x'^ — 3a1 — c^ + (?)xy'^ + \—a\{aia2X + 2c^)a;^ 

— aia2(— 3aj^ + 2c + c )a; — 2c (c — 2a;^)a; + c aia2j|/ 

+bib2x{a\x^ — c^)^ , (53c) 



2a/2( 



2 



{-aia26i(x^ - l)y^ 



■"^ (^2-^2)1/2(^2+ ^2)5/2 

+6^62(30^3;^ — 'ia^x + c'^a; + (?x + 2caia2)y*^ 
+6i[a^(3aia2a; + Sc^)x^ — aia2(4a^ — Tc"^ — 3c^)a;^ 

+2c(-4c^ai + c^ + 2aia2 + c^)^ - 0102 (&? + c^)\f 

-b2[a{x^ - Aa\{c^ - bl)x^ - 2c^aia2(2c2 - l)x^ 

+ {-Ac^al + c^ + 4c^ai + 3ai - c^)x + 2c^aia2]y^ 

— 61 [3af a2a;"^ + Aca\{a\ + c^)a;^ — aia2{'ia\ — 2(? — 7c^)x'^ 

—2c^(2a\ — c^ — (?)x — c^aia2]y 

+b2x{a'lx'^ — c'^){alx'^ + 2ca2aix + b\ + c^)^ . (53d) 

In the above expressions x and y are given by 



2x = ^P + {a + 1)2 + ^P + {a- 1)2, (54a) 

2y = ^P + {a + 1)2 - ^f2 + (tt _ 1)2^ (54b) 

where f = r/k and a = Zq/A;, with a > 1. 
When 62 = we have 

~ _ 16c^a2y2{[a(a;2 + 1) + ^^l + ^2)^2 _ ^2~2| 

™ ~ (x2-y2)[(aa; + c2)2 + 62^2]3 • (55) 

In order to analyze the behavior of -Dkn? is enough to consider the expression 

Do = [a{x^ + l)+x{l + c^)]^ -b^f^, (56) 
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that can be written as 



Do = a{l + c^)R+[a{a-l) + 2 + r]+a{l + c^)R4a{a + l) + 2 + r] 
+ ^R+R.[{c^ + If + a^f^ + a^ + 3)] + ^[(c^ + 1)^ + 5a'] 

+ U\' + 1 + a\f' + 2a' + 6)] + ]^a\^ + 1)^ + a\a' + 2)], 



(57) 



where R± = \/f'^ + (a ± 1)^. Since a(a =F 1) + 2 > for any a, from (157|) follows that D always is 
a positive quantity for Kerr-Newman fields and therefore the eigenvalues of the energy-momentum 
tensor are always real quantities. So we conclude that these disks can be interpreted, for all the values 
of parameters, as a matter distribution with currents and purely azimuthal pressure and without heat 

flow ng. 

In order to study the behavior of D when 62 7^ and of the other physical quantities associated 
with the disks, we shall perform a graphical analysis of them for charged and magnetized Kerr-NUT 
disks with h\ = 0.2, 62 = 0.9, charged and magnetized Taub-NUT disks with 62 = 0.9, and Kerr- 
Newman disks with bi = 0.2, for a = 2 and c = 1.0, 1.5, 2.0, 2.5, 3.0, as functions of f. In Fig. [TJ^a) 
we plot D for charged and magnetized Kerr-NUT disks and in Fig. [TJ^fo) we plot D for charged and 
magnetized Taub-NUT disks. We found that D takes negative values after certain value of f = fo. 
Therefore these disks always have heat flow beginning at the root fo- We also computed this function 
for other values of the parameters and, in all the cases, we found the same behavior. 

In Figs. [T](c) - lU^/i) we show the energy density e and the azimuthal pressure p^. We see that the 
energy density is in all cases always a positive quantity in agreement with the weak energy condition, 
whereas the pressure becomes negative for a value of f > fo in the case of charged and magnetized 
Kerr-NUT and Taub-NUT disks. We also see that the presence of electromagnetic field decreases the 
energy density at the central region of the disk and later increases it. The heat flow function q = kq 
is also represented in Figs, ^i) and Wij) for charged and magnetized Kerr-NUT and Taub-NUT 
disks. 

The electric charge density jj and the azimuthal current density j^, measured in the coordinates 
frame, are represented in Fig. [2l We observe that the electric charge density has a similar behavior 
to the energy density which is consistent with the fact that the mass is more concentrated in the 
disks center. We also computed this functions for other values of the parameters and, in all the cases, 
we found the same behavior. 

Now we will interpret the matter in the disks as two count err rotating streams of charged dust 
particles (CRM). First we consider the two streams of particles moving on electrogeodesics for charged 
and magnetized Kerr-NUT disks with bi = 0.2, 62 = 0.9, and Kerr-Newman disks with 61 = 2, for 
a = 1.2, and c = 1.0, 1.1, 1.3, 1.4, as functions of f. For charged and magnetized Taub-NUT fields 
we found a similar behavior to the charged and magnetized Kerr-NUT fields. 

In Figs, m^a) - ^d) we plot the angular velocities curves uj±. We see that these velocities are 
real and for charged and magnetized Kerr-NUT field 0;+ presents a nonsmooth behavior at certain 
value of f. In Figs. [3](e) - [3](/i) we have plotted the energy densities e±. We see that for charged 
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and magnetized Kerr-NUT disks e_ is everywhere positive whereas e+ becomes negative after certain 
value of f in violation of the weak energy condition. Therefore for charged and magnetized Kerr- 
NUT fields one finds that only the central region of the these disks presents a physically reasonable 
behavior. Thus these disks only model the inner portions of galaxies or acretion disks. However, 
since the surface energy density decreases rapidly one can to define a cut off radius at the point 
where the energy density changes of sign and, in principle, we consider these disks as finite so that 
we can use them to model the whole galaxy or acretion disk. By contrast for Kerr-Newman disks 
the counterrotating energy densities are always positive quantities. 

In Figs. [3](i) -EJ^O we plotted the electric charge densities a±. We also find that these quantities 
have a similar behavior to e±. In Figs. [3](?7i) -EJ^p) we have drawn the specific angular momenta h"^ 
for the same values of the parameters. We see that for charged and magnetized Kerr-NUT disks there 
is a strong change in the slope at certain values of f, which means that there is a strong instability 
there. We also find regions with negative slope where the CRM is also unstable. However, these 
disks models are stable in the central region and away of the center of the disks. Meanwhile for 
Kerr-Newman disks we find that the presence of electromagnetic field can also make unstable these 
disks against radial perturbations. Thus the CRM cannot apply for c = 3.0 (bottom curves). We 
also computed this functions for other values of 61 and 62 and, in all the cases, we found the same 
behavior. 

Second we consider the case when the two fluids move with equal and opposite tangential velocities 
(non-electrogeodesic motion) for charged and magnetized Kerr-NUT disks with a = 2, 61 = 62 = 0.1 
and c = 1.0, 1.5, 2.0, 2.5. In Fig. Hl^a) we plot the tangential velocity curves of the counterrotating 
streams, v^. We see that the quantity f ^ is also always less than the light velocity, but after of certain 
value of r becomes negative. Therefore these disks are also well behaved only in the central region of 
the disk. In Figs. \M,b) - \Myd) we have plotted the energy densities e± and the electric charge densities 
a±. We see that this quantities have a similar behavior to the previous cases. In Figs. W^e) andH)^/) 
we have drawn the specific angular momena h"^ and h'^ for the same values of the parameters. We 
find that these disks models are stable only in the central region of the disks. In the case of charged 
and magnetized Taub-NUT disks we also find that the physical quantities have a similar behavior to 
the charged and magnetized Kerr-NUT disks, and in the case of Kerr-Newman disks to the previous 
ones. 

Finally, the Figs. [T] - H] show that the two fluids are continuous in r which implies to have 
two particles in counterrotating movement in the same point in spacetime. So this model could be 
possible when the distance between streams (or between the counterrotating particles) were very 
small in comparing with the length r so that we can consider, in principle, the fluids continuous like 
is the case of counterrotating gas disks present in disk galaxies. 

5 Discussion 

We presented a detailed analysis of the energy-momentum tensor and the surface current density 
for electrovacuum stationary axially symmetric relativistic thin disks without radial stress and when 
there is heat flow. The surface energy-momentum tensor and the surface current density were ex- 
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pressed in terms of the comoving tetrad and explicit expressions were obtained for the kinematical 
and dynamical variables that characterize the disks. That is, we obtained expressions for the velocity 
vector of the disks, as well for the energy density, azimuthal pressure, electric charge density and 
azimuthal current density. 

We also presented in this paper the stationary generalization of the counterrotating model (CRM) 
for electrovacuum thin disks that had been previously only analyzed for the static case in [211 123] . 
We considered both counter rotation with equal and opposite velocities and counter rotation along 
electrogeodesics. A general constraint over the counterrotating tangential velocities was obtained, 
needed to cast the surface energy-momentum tensor of the disk in such a way that can be interpreted 
as the superposition of two counterrotating dust fluids. The constraint obtained is the generalization 
of the obtained for the vacuum case in [H], for disks without radial pressure and with heat flow, 
where we only consider counterrotating fluids circulating along geodesies. We also found that, in 
general, there is not possible to take the two counterrotating tangential velocities as equal and 
opposite neither take the two counterrotating fluids as circulating along geodesies. 

A family of models of counterrotating and rotating relativistic thin disks based on a charged and 
magnetized Kerr-NUT metric are constructed where we obtain some disks with a CRM well behaved. 
This solution has as limiting cases a charged and magnetized Taub-NUT solution and the well known 
Kerr-Newman solutions. For charged and magnetized Kerr-NUT and Taub-NUT fields we find that 
these disks always present heat flow, whereas for Kerr-Newman fields we shows that these disks can 
be interpreted, for all the values of parameters, as a matter distribution with currents and purely 
azimuthal pressure and without heat flow. Finally, the inclusion of radial pressure to these models 
is being considered. 
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Figure 1: The discriminant D (a) for charged and magnetized Kerr-NUT disks with bi = 0.2, 
62 = 0.9 and (b) for charged and magnetized Taub-NUT disks with 62 = 0.9, the energy density e (c) 
for charged and magnetized Kerr-NUT disks with bi = 0.2, 62 = 0.9, (d) for charged and magnetized 
Taub-NUT disks with 62 = 0.9 and (e) for Kerr-Newman disks with 61 = 0.2, the azimuthal pressure 
pip (/) for charged and magnetized Kerr-NUT disks with bi = 0.2, 62 = 0.9, (g) for charged and 
magnetized Taub-NUT disks with 62 = 0.9 and (h) for Kerr-Newman disks with 61 = 0.2, the heat 
flow function q = kq {i) for charged and magnetized Kerr-NUT disks with bi = 0.2, 62 = 0.9 and (j) 
for charged and magnetized Taub-NUT disks wi^h) 62 = 0.9, and a = 2.0, c = 1.0 (top curves), 1.5, 
2.0, 2.5, and 3.0 (bottom curves), as functions of f. 




r 
(a) 

Jip 



r 
(b) 

Jip 






r 

id) 



r 
if) 



Figure 2: The electric charge density j^ (a) for charged and magnetized Kerr-NUT disks with 
bi = 0.2, 62 = 0.9, (6) for charged and magnetized Taub-NUT disks with 62 = 0.9, and a = 2.0, 
c = 1.0 (axis f), 1.5 (top curves), 2.0, 2.5, and 3.0 (bottom curves), and (c) for Kerr-Newman disks 
with 61 = 0.2, and a = 2.0, c = 1.0 (axis f), 2.0 (top curve), 1.5, 2.5, and 3.0 (bottom curve), 
and the azimuthal current density j^ (d) for charged and magnetized Kerr-NUT disks with 61 = 0.2, 
62 = 0.9, (e) for charged and magnetized Taub-NUT disks with 62 = 0.9, and (/) for Kerr-Newman 
disks with bi = 0.2, and a = 2.0, c = 1.0 (axis f), 1.5 (top curves), 2.0, 2.5, and 3.0 (bottom curves), 
as functions of f. 
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Figure 3: We plot for electrogeodesic charged and magnetized Kerr-NUT fields with a = 1.2, bi = 0.2, 
62 = 0.9, and electrogeodesic Kerr-Newman fields with a = 1.2, 61 = 0.2, the angular velocities uj± 
(a), (6) for charged and magnetized Kerr-NUT disks and (c), (d) for Kerr-Newman disks, the energy 
densities e± (e), (/) for charged and magnetized Kerr-NUT disks and {g), {h) for Kerr-Newman 
disks, with c = 1.0 (top curves), 1.1, 1.3, 1.4 (bottom curves), the electric charge densities a± (i), (j) 
for charged and magnetized Kerr-NUT disks and (k), (/) for Kerr-Newman disks, with c = 1.0 (axis 
f), 1.1 (bottom curves), 1.3, 1.4 (top curves), th^|pecific angular momenta h'^ (m), (n) for charged 
and magnetized Kerr-NUT disks and (o), (p) for Kerr-Newman disks, with c = 1.0 (top curves), 1.1, 
1.3, 1.4, 3.0 (botton curves). 
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Figure 4: We plot for non-electrogeodesic charged and magnetized Kerr-NUT field with a = 2, 
hi =1)2 = 0.1, (a) the tangential velocity f^, (6) the energy densities e± for c = 1.0 (top curves), 1.5, 
2.0, 2.5 (bottom curves), (c), (d) the electric charge densities a± for c = 1.0 (axis f), 1.5 (bottom 
curves), 2.0, 2.5 (top curves), (e) and (/) the angular momentum specific h\ for c = 1.0 (top curves), 
1.5, 2.0, 2.5 (bottom curves), as functions off. 
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